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We study how the order of A*" independent random walks in one dimension evolves with time. 
Our focus is statistical properties of the inversion number m, defined as the number of pairs that 
are out of sort with respect to the initial configuration. In the steady-state, the distribution of the 
inversion number is Gaussian with the average (m) ~ N'^ /A and the standard deviation a ~ N^''^ /6. 
The survival probability, Sm{t), which measures the likelihood that the inversion number remains 
below m until time t, decays algebraically in the long-time limit, Sm ~ t~^"^ . Interestingly, there is a 
spectrum of N{N — 1) /2 distinct exponents Prn{N). We also find that the kinetics of first-passage in a 
circular cone provides a good approximation for these exponents. When is large, the first-passage 
exponents are a universal function of a single scaling variable, l3m{N) — )■ j3{z) with z — (m — {m})/a. 
In the cone approximation, the scaling function is a root of a transcendental equation involving 
the parabolic cylinder equation, D2p{—z) = 0, and surprisingly, numerical simulations show this 
prediction to be exact. 
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I. INTRODUCTION 

Consider the permutation 3142 of the four elements 
{1,2,3,4}. Three pairs: (1,3), (2,3), and (2,4) are in- 
verted in this permutation. The inversion number, de- 
fined as the total number of pairs that are out of sort, 
provides a natural measure for how "scrambled" a list 
of elements is. This basic combinatorial quantity 1-4] is 
helpful in many contexts. In computer science, the inver- 
sion number plays an important role in sorting and rank- 
ing algorithms [5] . Common on the web ( "customers who 
like . . . may also like . . ."), recommendations for books, 
songs, and movies use inversions to quantify how close 
the preferences of two customers are [a] ■ 

The number of inversions can also be used to measure 
how the order of a group of particles in one dimension 
changes with time. Figure [T] illustrates a space-time dia- 
gram of four diffusing particles. The number of inversions 
changes whenever two trajectories cross. Depending on 
the initial order of the two respective particles, a crossing 
may either introduce a new inversion or undo an existing 
one. Consequently, the inversion number either increases 
or decreases by one. Therefore, the inversion number 
equals the difference between the number of crossings of 
the first kind and the number of crossings of the second 
kind. 

Mixing dynamics has been extensively studied in the 
context of fluids 0, Q and granular materials Q, but 
much less attention has been given to mixing in the con- 
text of diffusion ^[il|. In this study, we consider an 
ensemble of N diffusing particles in one-dimension, a sys- 
tem that is widely used to model the transport of colloidal 
and biological particles in narrow channels [l^ [l^ . We 
use the inversion number to measure the degree to which 
particles mix. Clearly, a persistent small inversion num- 
ber indicates a poorly mixed system, while a large inver- 
sion number implies that the opposite is true. 

We first study how the distribution of the inversion 
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FIG. 1: Space-time diagram of a four-particle system. The 
circled +s and — s indicate whether the inversion number in- 
creases or decreases when two trajectories cross. Four out of 
the five crossings increase the inversion number, and accord- 
ingly, the inversion number increases from m = to m = 3. 



number evolves with time. We find that there is a 
transient regime in which the average inversion num- 
ber as well as the standard deviation in this quantity 
both grow as the square-root of time. The distribu- 
tion of the inversion number is stationary beyond this 
transient regime. When the number of particles is suffi- 
ciently large, the probability distribution function is al- 
ways Gaussian, whether in the transient regime or in the 
steady-state. 

Our main focus is first-passage properties [l3| of the 
inversion number. We ask: what is Sm,{t), the probabil- 
ity that the inversion number remains smaller than m up 
to time t. For small values of m, the survival probabil- 
ity Sm measures the likelihood that the particles remain 
poorly mixed throughout the evolution. Generally, the 
probabilities Sm decay as a power law at large times, 
Sm ^ t"^'" . In general, there is a broad spectrum of 
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N{N -l)/2 distinct exponents, /32, • ■ ■ , 13n{n-i)/2}, 
that governs the asymptotic decay of the survival prob- 
abilities. 

We heavily use first-passage kinetics of a single parti- 
cle that diffuses inside a circular cone [TEl - fr^ to under- 
stand the asymptotic behavior of Sm- We first utilize 
two-dimensional cones to obtain the first-passage expo- 
nents for a three-particle system exactly. Furthermore, 
we employ circular cones in A'^ — f dimensions and find 
good approximate values for the first-passage exponents. 

The cone approximation correctly predicts that when 
the number of particles is large, the exponents become a 
universal function, Pm{N) — >■ f}{z), of the scaling variable 
z = [m — {m))/(T. Here, (m) and a are the average 
and standard deviation of the distribution of inversion 
number, at the steady-state. Interestingly, our numerical 
simulations show that the cone approximation yields the 
exact scaling function as a root of a transcendental 
equation involving the parabolic cylinder function. 

The rest of this paper is organized as follows. In section 
II, we introduce our basic system and define the inversion 
number. Stationary and transient properties of the distri- 
bution of the inversion number are discussed in sections 
III and IV, respectively. In section V, we use the cone 
approximation to understand first-passage properties of 
the inversion number. Scaling and extremal properties 
of the first-passage exponents are the focus of section VI. 
We conclude in section VII. 



II. THE INVERSION NUMBER 

Our goal is to characterize how the order of an ensem- 
ble of diffusing particles changes with time. We conve- 
niently use an ordinary random walk [l8l - [2^ to model 
the trajectory of a diffusing particle [2l|. Our system in- 
cludes N identical particles that move on an unbounded 
one-dimensional lattice. The particles are completely in- 
dependent: at each step one particle is selected at ran- 
dom and it moves, with an equal probability, either to 
the left, a; — >■ a; — I, or to the right, x — > a; + I. After 
each elementary step, time is augmented by the inverse 
number of particles, t —i' t + 1/N, so that each particle 
moves once per unit time. 

We index the particles according to their initial po- 
sition with the leftmost particle labeled n = 1 and the 
rightmost particle labeled n = N (Figure[T]). Let Xn{t) be 
the position of the nth particle at time t. By definition, 



a;i(0) < X2(0) < • • • < a;Ar_i(0) < xn{0), 



(1) 



but the initial order unravels with time. Consider, for 
example, the four-particle system illustrated in Figure 
[TJ The particles reach a state where X3 < X2 < xi < X4 
with three pairs, (I, 2), (I, 3), and (2, 3) being out of sort 
compared with time t = 0. In general, a pair of particles 
for which Xi (t) > Xj (t) and i < j constitutes an inversion. 



Formally, the total number of inversions, m, is given by 

N N 

^w = E E Q{^^{t)-^At))■ (2) 

Here, Q{x) is the Heaviside step function; 9(x) ~ I for 
X > and Q{x) = for x < 0. The total number of 
pairs is Af = {2) > hence, the variable m is within 
the bounds < to < M with 



M = 



N{N - 1) 



(3) 



The inversion number is minimal, m — 0, when the order 
is exactly the same as in the initial configuration, and it 
is maximal, to = M, when the order is the mirror image 
of the initial state. 

The inversion number changes whenever two trajec- 
tories cross (Figure [T]). A crossing either adds a new 
inversion or removes an existing one. Thus, we may as- 
sign a positive or a negative "charge" to each crossing as 
illustrated in Figure [TJ The inversion number, m{t), is 
simply the sum of all of the charges up to time t. 



III. THE MAHONIAN DISTRIBUTION 

We first discuss basic statistical characteristics of the 
inversion number including the average, the variance, and 
more generally, the probability distribution function. At 
large time i, each random walk explores a region of size 
^/i, and the probability of finding the particle at any lat- 
tice site inside this region is effectively uniform. This 
simple fact already implies that memory of the initial 
position fades with time. We thus expect that after suffi- 
cient time elapses, there is no memory of the initial order, 
and the order of the particles is completely random. 

To understand statistics of the inversion number for 
randomly ordered particles we consider the set of all A^! 
permutations of the N elements {I, 2, . . . , iV}. In the 
random state, each permutation of these elements occurs 
with probability 1/iV!. The probability Pm{N) that the 
inversion number equals to for a random permutation is 
well known as the Mahonian distribution in probability 
theory [H, 01 • We highlight key features of this proba- 
bility distribution as it plays a central role in our study. 

Let Qm{N) = NlP,n{N) be the number of permuta- 
tions of N elements with exactly to inversions. For exam- 
ple, when TV = 3, one permutation (123) is free of inver- 
sions, there are two permutations with one inversion (213, 
132), two permutations with two inversions (312, 231), 
and a single permutation with three inversions (321). 
Hence, Qo(3) = ^3(3) = 1 while Qi(3) = Q2(3) = 2. We 
list the distribution of the inversion number for < 4, 
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•(1) iV=l, 

(1,1) N^2, 

(1,2,2,1) iV-3, 

(1,3,5,6,5,3,1) iV = 4. 
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Of course, Pm{N) is nonzero if and only if < m < M. 

Since the mirror image of a configuration with m inver- 
sions necessarily has M — m inversions, the probability 
distribution satisfies Pm = PM-m- Hence, the distri- 
bution is symmetric about m = M /2, and the average 
(to) = '^P^n is simply 



(to) 



N{N -1) 



(4) 



Therefore, the average grows quadratically with the total 
number of particles when N ^ 1. 

The Mahonian distribution satisfies the simple recur- 
sion relation 



(5) 



1=0 



with /'m(l) — Sm.o- This recursion reflects that every 
permutation of iV elements can be generated from a per- 
mutation of A'^ — 1 elements by inserting the iVth element 
in any of the N possible positions. Depending on where 
this last element is added, the number of inversions in- 
creases by an amount Am = 0, 1, . . . , TV — 1. 
Let us now introduce the generating function. 



M 



(6) 



For instance, T(.s, 1) = 1, T(s,2) = (1 -f s)/2! and 
J'(s, 3) = (l + s)(l + s + s^)/3!. In general, the generating 
function is given by the product [5| 



1 ^ 



+ s + + 



+ s 



n-lN 



(7) 



as also follows from the recursion ([5]). 

We can confirm the average (jlj by differentiating the 
generating function once, CP'(s = 1) = (m), where 
the prime represents differentiation with respect s. By 
differentiating the generating function twice and using 
J"'(s = 1) = {m{m — 1)), we obtain the variance Q, 



(to^) 



N{N -1){2N + 5) 
72 ' 



(8) 



This expression is obtained from = TX"' There- 

fore, the standard deviation is rather large, a ~ 
when iV > 1. 

The mean ([4]) and the variance ([8]) fully specify the 
probability distribution function for an asymptotically 
large number of particles. The Mahonian distribution 
becomes a function of a single variable, Pm{N) — >■ ^{z), 
with the scaling variable 



(9) 



The probability distribution function, $(2;), is normal, 
that is, a Gaussian with zero mean and unit variance 

am, 



$(z) 



1 



exp 



(10) 



To see that the central limit theorem applies, we con- 
vert the generating function into a Fourier transform, 
and then show that the Fourier transform is Gaussian in 
the large- TV hmit [HI. 

The variable z is a more transparent measure in the 
following sense. A value of z of order one implies fairly 
random order. Indeed, according to the normal distribu- 
tion (jlOp . the inversion number falls within three stan- 
dard deviations from the mean, \z\ < 3, with probability 
0.997. A large value, \z\ ^ 1, indicates that the parti- 
cle order strongly resembles the initial configuration (if 
z > 0) or its mirror image (if z < 0). 

IV. TRANSIENT BEHAVIOR 

By definition, the inversion number is zero initially, 
to(0) = 0. At least partially, the initial order is pre- 
served in the early stages of evolution, and the number 
of inversions must be substantially lower than Q . 

We consider the natural initial condition where the 
particles occupy N consecutive lattice sites: 2:^(0) = i, 
for all i — 1,2, . . . , N. Early on, particles "interact" 
only within their local neighborhood. The interaction 
length, i, grows diffusively with time, £ ~ \/t. On 
this length scale, particles are well-mixed, and according 
to ([U, the number of inversions per particle is propor- 
tional to the number of interacting particles, i. Hence, 
the average number of inversions grows according to 
(TO a consequence. 



{m{t)) 
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iV^ < t. 



(11) 



when ^ 1. The two expressions match at i ~ N'^, 
a diffusive time scale that can be viewed as the mixing 
time. Therefore, there is a transient regime in which 
the inversion number grows as the square-root of time, 
followed by a steady-state, in which the average is given 
by Q. 

We obtain the variance using a similar heuristic ar- 
gument. According to ([H|), the variance per parti- 
cle is quadratic in the number of interacting particles, 
~ Nl'^ . Therefore, cr^ ~ TV i in the transient regime. 
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{const X 
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1< t < A^2^ 

< t. 



(12) 



As expected, the transient behavior matches the steady- 
state behavior (jS]) at the diffusive time scale t ~ iV^. 
Like the average, the standard deviation also grows as 
the square root of time. 
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FIG. 2: The normalized average {m{t))/[N{N - l)/4] versus 
time t. The results correspond to an average over 10^ inde- 
pendent realizations of a system with A'^ = 10'' random walks. 
Also shown for reference is a line with slope 1/2. 



As shown in Figure [21 results of numerical simulations 
confirm the scaling behavior (jlll) . Moreover, the numer- 
ically measured average matches the steady-state value 
corresponding to the Mahonian distribution. We also 
verified that the stationary distribution is Gaussian with 
the variance dS]). 

The simulations also show that the time-dependent 
distribution of inversion number, pm{N,t), is Gaussian 
throughout the transient regime (Figure E]): 



1 



exp 



(to — (m(t)))^ 



2<7^{t) 



(13) 



This behavior provides further support for our heuristic 
argument. Indeed, if the particles are well-mixed locally, 
then the distribution of the number of inversions per par- 
ticle is Gaussian, and as the sum of N Gaussian variables, 
the total inversion number must also have a Gaussian dis- 
tribution. 

We used two different algorithms to simulate the diffu- 
sion process. In the naive algorithm, we randomly select 
a particle and move it to a randomly-chosen neighboring 
site. We increase time by l/iV after each jump. To cal- 
culate the inversion number, we use the formula but 
since this enumeration requires 0{N'^) operations, this 
simulation method is inefficient at large N. 

To overcome this difficulty, we introduced a variant 
where each lattice site may be occupied by at most one 
particle. At each step we pick one particle at random 
and attempt to move it by one site. This move is always 
accepted if the neighboring site is vacant, but otherwise, 
it is accepted with probability 1/2. In the latter case, 
we merely exchange the identities of the respective parti- 
cles, and as appropriate, update the inversion number by 
either -1-1 or —1. In our implementation, there are two 
arrays: the first lists the particle positions, in order, and 
the second lists the original position of each particle in 
the first list. This algorithm has a fixed computational 
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FIG. 3: The distribution of the inversion number in the inter- 
mediate time regime 1 <C t A'^^. Shown is the distribution 
Pm = Pm{N,t) versus the variable (m — {m{t))^ /cr{t). The 
results are from 10^ independent realizations of a system with 
A*' — 10^ random walks. The distribution is shown at times 
t — 10^ (diamonds), t = 10^ (squares), and t = 10* (circles). 
Also shown for reference is the normal distribution H13[). 



cost per step, and it automatically keeps track of the 
inversion number. We rely on the fact that in one di- 
mension, noninteracting random walks are equivalent to 
random walks that interact by exclusion [2J-|27|. Still, 
we verified that the two algorithms yield essentially the 
same results. We utilized the naive algorithm to simulate 
small systems with iV < 10, but otherwise, we used the 
efficient algorithm. 



V. FIRST-PASSAGE KINETICS 

We have seen that the inversion number, which grows 
quadratically with the number of particles, can be quite 
large. Yet, if the mixing is poor and the particle trajecto- 
ries rarely cross, the inversion number remains small. To 
quantify how common such a scenario is, we study first- 
passage kinetics [H, [28|. In particular, we ask: what 
is the probability, Sm{t), that the inversion number re- 
mains smaller than to until time t. This "survival" prob- 
ability is closely related to the first-passage probability as 
[—dSm/dt] X dt is the probability that the inversion num- 
ber reaches to for the first time during the infinitesimal 
time interval (i, t + dt). 

The quantity Si is the probability that the original 
order is perfectly maintained, or equivalently, the likeli- 
hood that none of the trajectories cross. This survival 
probability decays as a power law, with a rather large 
exponent, 

^1 ^ t-^(^-^)/\ (14) 

in the long-time limit (29l - [33 | . Our goal is to understand 
how this asymptotic behavior changes as the threshold 
TO increases. 



When N = 2, the separation between the two random 
walks itself undergoes a random walk. Hence, Si is equiv- 
alent to the survival probability of a one-dimensional ran- 
dom walk in the vicinity of a trap, and Si ^ t^^/^ in 
agreement with ([T4|) . 

When iV = 3, we conveniently map the three random 
walks onto a single "compound" random walk in three 
dimensions with the coordinates {xi,X2,X3). To find 5*1, 
we require that the compound walk remains inside the 
region xi < X2 < x^. We may view the boundary of this 
region as absorbing, and then, (t) equals the likelihood 
that that the compound walk survives at time t. The ab- 
sorbing boundary forms a wedge because it is defined by 
the intersection of two planes, xi = X2 and xi ~ x^- Gen- 
erally, the survival probability of a particle that diffuses 
inside an absorbing wedge decays algebraically. 



(15) 



where V = a/ir is the normalized opening angle |35| . 
(The opening angle Q < a < n is the angle between 
the wedge axis and the wedge boundary.) Alternatively, 
< y < 1 is the fraction of the total solid angle en- 
closed by the wedge. The region xi < X2 < x^ occupies 
a fraction Vi = ^ = g o f space and hence, 5*1 ~ 
as also follows from (|14l) . To find 5*2 and 5*3, we note 
that the regions in which the compound walk is allowed 
to move are always wedges (the three planes xi = X2, 
xi = X3, and X2 = X3 divide space into six equal wedges 
[3^.) Moreover, the fraction of total solid angle enclosed 
by the absorbing boundaries is given by the cumulative 
distribution of inversion number: 1^2 = + fr = 5 ^-nd 
V3 = ^-|-|^-|-|i=|. Hence, all three survival proba- 
bilities decay algebraically with time [l^ , 
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(16) 



and there are three distinct first-passage exponents. 

The asymptotic behaviors (|16p suggest that all of the 
survival probabilities decay algebraically. 
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(17) 



in the long-time limit. Moreover, there is a large family 
of exponents 



Pi > 132 > ■ • • > I^N(N~l)/2i 



(18) 



that characterizes the power- law decay ([T7)) . We stress 
that the exponents depend on two variables, the thresh- 
old m and the number of particles iV, /3 = /?,„(iV). We 
already know the exact values /3i(3) = 3/2, /32(3) = 1/2, 
and /33(3) = 3/10 as weh as ^i(7V) = N{N 

Our numerical simulations confirm that indeed, there 
is a large spectrum of exponents. As shown in Figure 31 
there are six decay exponents when = 4. Table I lists 
the numerically measured values /3„i, obtained from the 
local slope din Sm I dint. 

In general, the compound walk is confined to a certain 
"allowed" region of space. This region is bounded by mul- 
tiple intersecting planes of the type Xi = Xj with i j, 




10 1234567 

10 10 10 10 10 10 10 10 

t 



FIG. 4: The survival probability Smit) versus t for a four- 
particle system. Shown are the quantities S2 (bottom curve), 
S3, . . ., Se (top curve). The number of independent Monte 
Carlo runs varies from lO" for the slowest decay to 10^^ for 
the fastest decay. 



and generally, this unbounded domain has a complicated 
geometry. The boundary of this region encloses a frac- 
tion Vrfi{N) of the total solid angle. On combinatorial 
grounds alone, we conveniently deduce that this fraction 
is given by the cumulative Mahonian distribution 



(19) 



1=0 



Since the Mahonian distribution is symmetric, we have 
V,n + Vaz+i-to = 1- To evaluate VmiN), we expand the 
generating function ([7]), and for to < 4, we have 5] 



1 

N 



m — 1 
2 



^{N -1){N + 2) m = 3, 

1{N + 1){N^ + 2N - 6) m = \. 



' (20) 



We have seen that the allowed region is a wedge when 
A^ = 3. To obtain an approximation for the first-passage 
exponents, we follow an approach that proved useful in 
other first-passage problems involving multiple random 
walks and replace the boundary of the allowed region 
with a suitably chosen cone in A^ — 1 dimensions [s^l ■ An 
unbounded cone with opening angle a occupies a fraction 
VioL) of the total solid angle, given by 
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/;d0(sin0) 
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Ar-3 



(21) 



In d dimensions, we have dil c>c (sva^OY-^^dB where 
is the solid angle and is the polar angle in spherical 
coordinates. In the cone approximation, we require 
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2.67100 
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0.64975 


0.39047 


0.24517 


0.14988 
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1.39 


0.839 


0.455 


0.275 


0.160 



TABLE I: The six first-passage exponents for a four-particle 
system. The values /J™, are from the Monte Carlo simulation 
results shown in Figure H The values ^^""^ were obtained us- 
ing the cone approximation, specified in Eqs. (|19|l - (|23[ ). The 
cumulative Mahonian distribution, Vm, and the opening an- 
gle, am, are listed as well. 



with Vjn given in (1191) . 

In a cone, the first-passage exponent (3 = j3{a) de- 
creases as the opening angle a increases. In particular, 
/? = Ti/Aa in two dimensions, and /? = (tt — a) /2a in four 
dimensions. Generally, however, /3 is the smallest root of 
the following transcendental equation involving the asso- 
ciated Legendre functions [ssll of degree 2/3-1-7 and order 



Poa , .^(cosa) 



2/3-1-7 I 



N odd, 
N even. 



(23) 



Regardless of the dimension, the surface of a cone with 
a = tt/2 is a plane, and hence, /3(7r/2) = 1/2. 

For example, to find /3i (4) , we first determine the frac- 
tion ^1(4) = jj- = ^ using ([T9l) . Then, we calculate 
the opening angle a = 0.41113 using equations (PT|) - (P^ 
and finally determine the exponent /3i(4) = 2.67100 as 
the appropriate root of equation ([23| . By construction, 
the cone approximation is exact for three particles. This 
approach gives a useful approximation to the six first- 
passage exponents when = 4 (Table I). Remarkably, 
the cone approximation continues to be a good approxi- 
mation as the number of particles increases (Figure [5]) . 



VI. THE SCALING FUNCTION 

We are especially interested in the behavior when the 
number of particles is large. Let us first evaluate the 
cumulative Mahonian distribution in the large- A'^ limit. 
Since the Mahonian distribution is normal, the cumula- 
tive distribution is given by the error function, 



VUN) 



(24) 



when N ^ 00. Here, z is the scaling variable defined 

jp 

into 



in ^ and erf(^) = (2/^/7:) cxp{-u'^)du. To obtain 
Eq. (1211), we substitute ^ into ^ and convert the 
sum into an integral. Equation (|24p is relevant in the 
limit iV — >■ (X), m — >■ (X) with the scaling variable z finite. 

Next, we evaluate the solid angle enclosed by an un- 
bounded cone when the dimension is large. The domi- 
nant contribution to the integral in (|2ip comes from a 
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FIG. 5: The first-passage exponent versus m for A'^ = 
4, 5, 6, 7. Shown are simulation results (circles) and the out- 
come of the cone approximation (squares). 



narrow region of order 1 /v N centered on a = 7r/2 where 
the integrand is Gaussian, 



(sin( 



\N~2 



-N{Ti/2-ey/2 



Using exp [-N{Tr/2 - Of /2\de ^27r/A^, we find 
that the fraction V{a) has the scaling form 



V{a,N) 



f(-y 



(25) 



with y — {cos a)\'N. In writing this equation, we used 
the facts that cosa~7r/2 — a and erf(^) = — crf(— ^). 
Equation (|25|) holds in the limit A^-^oo, a— >7r/2, with 
the scaling variable y finite. 

Asymptotic analysis of equation shows that the 
exponent (3 {a) adheres to the scaling form [l3| 



(26) 



Pia,N) l3{y) with y = (cos a)VN 



in the limit N ^ 00, a ^ tt /2 with the scaling variable 
y finite. The scaling function, f5{y), is specified by the 
transcendental equation D2p{y) — 0, where D^, is the 
parabolic cylinder function of order v The smallest 
root is the appropriate one (l7| . 

By comparing equations (^^ and (PSI) . we find our 
main result: the first-passage exponent depends on a sin- 
gle scaling variable. 



p^[N)^P{z) with z 



(27) 



in the large- A^ limit. We reiterate that the average (m) 
and the standard-deviation a correspond to the steady- 
state values (01) and dHl), respectively. Using y = —z, the 
scaling function f3{z) is the smallest root of the transcen- 
dental equation 
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FIG. 6: The exponent /3 versus the scahng variable z, shown 
using: (a) a linear- linear plot and (b) a linear- log plot. The 
simulation results are from Monte Carlo runs with A*" = 50 
(diamonds), = 100 (squares), and A'^ = 200 (circles) parti- 
cles. The number of independent realizations varies from 10* 
for slow first-passage processes to 10* for fast one. The solid 
line shows the theoretical prediction 



involving the parabolic cylinder function. When /3 is a 
half-integer, the parabolic cylinder function is related to 
the Hermite polynomials |38| and using this equivalence, 
we have /3(0) = 1/2, /3(-l) = 1, and /3(-\/3) 3/2. 

Our numerical simulations (Figure 15]) confirm that the 
exponents /3m{N) have the scaling form (f27|) . Interest- 
ingly, the simulations strongly suggest that the scaling 
function predicted by the cone approximation is exact. 
We note that the convergence to the infinite-particle limit 
is very fast for positive z, but much slower for negative 

With the power-law decay (|17p , the mean first- passage 
time diverges whenever /? < 1, but it is finite otherwise. 
Since /3(z = —1) = 1, the time required for the inversion 
number to reach one standard deviation from the mean 
is infinite, on average. Regardless of the threshold z, 
there is a considerable chance that the random walks 
are poorly mixed because the survival probability decays 
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6 
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Pi 


3 
2 


3 


5 


15 
2 


21 
2 


14 


^conc 


3 
2 


2.67100 


4.08529 


5.73796 


7.62336 


9.73686 


^conc 

Pm 


3 

10 


0.14988 


0.061195 


0.019895 


0.0050713 


0.0010266 



TABLE II: The largest exponent, ^J""^ and the smallest 
exponent, Pm'^", obtained using the cone approximation for 
N < 8. Also listed for reference, is the exact value /3i. 



algebraically. 

The scaling behavior is remarkable for a number 
of reasons. First, the form of the scaling variable, 
z = {m- N^/4:)/{N^^'^/6), is quite unusual. Second, 
there are roughly first-passage exponents and nu- 

merical evaluation of this large spectrum is daunting. 
Yet, the scaling form ((27|) gives the range of parame- 
ters for which /3 is of order one, and hence, numerically 
measurable. (It is difficult to measure a vanishing expo- 
nent, /3 -> 0, or a divergent exponent, f3 oo.) Last, 
the emergence of scaling laws for a family of scaling expo- 
nents is also intriguing. Typically, in Statistical Physics, 
the opposite is true as one or two scaling exponents char- 
acterize a scaling law (soj . 

The extremal behaviors of the roots of the transcen- 
dental equation ([^5)) are derived in ref. (iTj . 



yiVS^exp (-^V2) 



— oo, 

CO. 



(29) 



The first-passage exponent is algebraically large if z is 
large and negative, but it is exponentially small if z is 
large and positive. 

The exponential decay in (I29p implies that it is ex- 
tremely unlikely that the initial order is perfectly re- 
versed. The smallest exponent /3m characterizes the 
probability Sm that the order of the walkers does not 
turn into the mirror image of the initial state, that is, the 
probability that the compound walk remains in the exte- 
rior oi the so-called "Weyl chamber" xi < X2 < ■ ■ ■ < 
[H-lll. This domain has Vm = 1-^ , and Table II lists 
the outcome of the cone approximation for small N. To 
find the outcome of the cone approximation at large N, 
we first estimate the opening angle, tt — a ~ e/N by using 
Eq. dSI]) and the Stirhng formula iV! ~ V^^{N/e)^. 
From the asymptotic behavior for wide cones at large di- 
mensions, /3 ~ yjN/%-K{TT - a)^~^ [l3l, we conclude ^ 



M 



(30) 



This value is extremely small, decaying roughly as the 
inverse of a factorial, and it is impossible to measure 
such a minuscule quantity using numerical simulations. 

The largest exponent describes the probability that the 
particles maintain the initial order or that the compound 
walk remains in the interior of the Weyl chamber with 
Vi — ^ . Table II compares the outcome of the cone ap- 
proximation with the exact value /3i = N{N — l)/4. The 
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quality of the cone approximation worsens as N grows. 
Nevertheless, the cone approximation is qualitatively cor- 
rect. By substituting the opening angle a c:^ e/N into 
the thin-cone asymptotic behavior /3{a) ~ 7Va~^/4 [17], 
we find Pi ~ iV^/(4e). This expression captures the 
quadratic growth of the exponent. Remarkably, the cone 
approximation is exact inside the scaling region, but it is 
only approximate outside this region. 

VII. CONCLUSIONS 

In summary, we used the number of pair inversions 
to measure the one-dimensional mixing of independent 
diffusing trajectories. A high inversion number typifies 
strong mixing whereas a persistent small inversion num- 
ber indicates poor mixing. In the steady-state, the dis- 
tribution of inversion number is given by the well-known 
Mahonian distribution, and consequently, it is Gaussian 
when the number of particles is large. Preceding the 
steady-state is a transient regime in which both the av- 
erage inversion number and the standard deviation grow 
diffusively with time. 

We focused on first-passage statistics and showed that 
the probability that the inversion number does not ex- 
ceed a certain threshold decays as power law with time. 
Moreover, we found that a large spectrum of decay expo- 
nents characterizes the asymptotic behavior. When the 
number of particles is large, the exponents obey a uni- 
versal scaling function. The scaling variable equals the 
distance between the threshold inversion number and the 
average inversion number, measured in terms of the stan- 



dard deviation. 

The cone approximation, which replaces the region in 
which the compound random walk is allowed to move 
with an unbounded circular cone, plays a central role in 
our analysis. This approach is exact for three particles, it 
produces very good estimates in higher dimensions, and 
remarkably, this framework yields the exact scaling func- 
tion. The cone approximation gives lower bounds for the 
decay exponents because, among all unbounded domains 
with the same solid angle, the perfectly circular cone 
maximizes the survival probability 37i, .4&. AL] ■ The cone 
approximation is useful in answering other first-passage 
questions such as the probability that the nth rightmost 
random walk does not cross the origin and the probability 
that the original rightmost particle always remains ahead 
of at least n other particles [s^j. In both cases, there are 
as many exponents as there are particles, and curiously, 
the circular cone framework produces the scaling func- 
tion governing the first-passage exponents approximately 
in the first case and exactly in the second case. 

Understanding when the cone approximation is exact 
and when it is approximate is an interestiiig challenge, 
with implications well beyond first-passage |42h44| . The 
first-passage exponent is directly related to the lowest 
eigenvalue of the Laplace operator, and therefore, we con- 
clude that the lowest eigenvalue of the Laplacian similarly 
obeys scaling laws in high dimensions. The shape of the 
scaling function depends on the actual geometry [isj . 

I thank Paul Krapivsky and Timothy Wallstrom for 
useful discussions. This research is supported by DOE 
grant DE-AC52-06NA25396. 
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